Fluctuation Theorem in a Quantum-Dot Aharonov-Bohm Interferometer 
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In the present study, we investigate the full counting statistics in a two-terminal Aharonov-Bohm 
interferometer embedded with an interacting quantum dot. We introduce a novel saddle-point solu- 
tion for a cumulant-generating function, which satisfies the fluctuation theorem and accounts for the 
interaction in the mean-field level approximation. Nonlinear transport coefficients satisfy universal 
relations imposed by microscopic reversibility, though the scattering matrix itself is not reversible. 
The skewness can be finite even in equilibrium, owing to the interaction and is proportional to the 
asymmetric component of nonlinear conductance. 
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Microscopic reversibility is a key ingredient in deriv- 
ing the Onsager relation and has played a fundamental 
role in establishing the linear response theory Re- 
cently, microscopic reversibility was used to develop a 
new relationship that would be valid beyond the linear 
response regime. This is now known as the fluctuation 
theorem (FT) |2j. The FT relates probabilities between 
positive and negative entropy productions; provides a 
precise statement for the second law of thermodynamics; 
and remarkably, reproduces the linear response theory, 
the Kubo formula and the Onsager relation Q. 

In the last few years, full counting statistics (FCS) 
has been recognized as a suitable framework for the FT 
in quantum transport [!, 0, H, @, 0, 0. FCS provides 
a comprehensive statistical properties for charge trans- 
port far from equilibrium 0, [l(| [H|. It addresses the 
probability distribution P(q) of the charge q which is 
transmitted during time t, and its cumulant generat- 
ing function (CGF): = lim T ^oo lnZ(A)/r, where 
Z{\) = J2 q P(q) e igA and A is called the counting field 0. 
Recently, the FT was generalized to the quantum trans- 
port regime in the presence of interaction and a magnetic 
field B [5|. For two-terminal systems, the FT is 



JF(A; B) = F{-\ + iA;-B). 
P(q;B) = P(-q;-B)ei A , 



(1) 
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where A is the affinity A — V/T, the ratio between voltage 
V and temperature T (e = h= ks = 1)- One important 
consequence from ([T]) is the universal relations among 
transport coefficients 0. The transport coefficient L is 
obtained by expanding the current cumulant with respect 
to A: 
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where I = q/r. The FT Q leads to the Kubo for- 
mula L\ = Lq/2 and the Onsager relation L\ _ = 



0, where L™ 1± = ± L^(-B) is the sym- 

metrized/antisymmetrized transport coefficient. Fur- 
thermore, nontrivial relations among higher-order coef- 
ficients are obtained 0: 

L\ ^= -L\_= -Lq _ , L\ + = l\ + , Z/q + = 0. (3) 

This is significant in that the skewness Lq _ can be finite 
even in equilibrium and proportional to the asymmetric 
component of nonlinear conductance L\ _ as well as the 
linear response of noise L\ _ . 

However, there are some controversies regarding the 
validity of Eq. Q} Q. When interacting mcsoscopic con- 
ductors possess no mirror symmetry, the nonequilibrium 
charge accumulation inside the conductor is not sym- 
metric in the magnetic field Then, the potential 
landscape generated by the nonequilibrium charge ac- 
cumulation is not symmetric either. This implies that 
the S'-matrix is not reversible with respect to a magnetic 
field Slr(B) ^ Srl(-B), which generates the magnetic 
field asymmetric component of nonlinear conductance 
and 'violates' the Onsager relation 0, [H, O, (3, Ell- 
in this respect, the FT |T]) is counterintuitive. In fact, 
in Ref. , the Mach-Zehnder interferometer is suggested 
as a counter example of the FT Hence, it is neces- 
sary to give examples preserving the FT |T|) . To this end, 
we consider a two-terminal Aharonov-Bohm (AB) inter- 
ferometer embedded with a quantum dot (QD) [inset in 
Fig. H](a)] |16j| . We introduce a novel saddle-point solu- 
tion of CGF, which realizes the FT ([1]) and the lack of re- 
versibility in the S'-matrix simultaneously. It is achieved 
by introducing the 'counting field of the dot charge' in 
addition to the dot potential, which are functions of A. 
The solution accounts for nonequilibrium charge accu- 
mulation and current fluctuations in the Hartree-level 
approximation. We will also calculate Eq. explicitly 
and show that the equilibrium skewness is a consequence 
of Coulomb interaction. 

CGF of QD AB interferometer - The system consists 
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of left (L) and right (R) leads, two arms and a QD. Elec- 
trons can travel through the QD and the lower reference 
arm [inset in Fig. [T](a)]. The total Hamiltonian is 



H 



E 

r=L,R 



H r + H D + H T + H r 



ef, 



(4) 



where the on-site Coulomb interaction U in the QD is 
accounted for by Hjj = i e od\d a + Ud^d^d^dy. 

The operator d a annihilates an electron with spin a. 
The leads are modeled by H r = J2kcr £ rka a\. kcT a r k<j, 
where a r ka annihilates electrons in the lead r with spin 
a and wave vector k. The tunneling and the reference 
arm are described as ifp = J2 r ka *r d\,a r k a + H.c. and 
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H.c. . The magnetic 



field B pierces through the ring, and the electrons ac- 
quire the AB phase cf>, which satisfies <j)(B) = —<fi(—B). 
The initial density matrices at both leads are assumed to 
have an equilibrium distribution with the chemical po- 
tential fiL = V/2 and //r = —V/2. 
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FIG. 1: (a) Aharonov-Bohm (AB) phase dependent nonequi- 
librium charge accumulation, (inset) Quantum-dot AB Inter- 
ferometer. Mirror symmetry along the horizontal axis (dotted 
line) is absent, (b) Demonstration of the fluctuation theo- 
rem (2). The inset shows probability distributions for posi- 
tive and negative magnetic fields. Parameters: FL/r = 0.25, 
r fl /r = 0.75, trcf = 0.25, e D =0, and U=T=V=Y. 

We confine ourselves to high temperature and treat 
the interaction in a mean-field level approximation [l9| . 
In order to find a proper saddle-point solution out of 
equilibrium, we employ the real-time path integral ap- 
proach [l7j |. The characteristic function, which is the 
partition function in the Keldysh formalism, reads Z = 



J V[a* ka , d* , a rka , d a ] exp (i J c dtC(t)), where C is the 
closed time-path. The Lagrangian is given by 

£ = £ a *ka( id t - s r ka) <Vfco-+y^<C ( id t -e D )d a 

rka <J 

— ^^(t r e lipr d^cirka + c.c.) — U d*d-[d^di 
£ (t LR e-^^a* Rka a L k, + c.c.) , 

kk' a 

where the phase on the upper/lower branch of the 
closed time-path tp r ± is related to the counting field 
as (p r ± — ±X r /2. We introduce the auxiliary dot- 
potential v„ via the Stratonovich-Hubbard transforma- 
tion: Ud*d]-d*di — ► J2a v a d%d a — v^vJU [II]: 



5 = I V[v a ] Z (v a ) exp ( - / dtv^t) Vl (t) 



where Zq(vo) is the Keldysh partition function for the 
noninteracting case U = with a shift in the QD level 
erj — * e_D+i> CT (t) f° r spin a. Although we limit ourselves 
to the time-independent stationary solution in the non- 
magnetic phase [ISfl , we allow different dot-potentials for 
upper and lower branches of C [2fj| : 



v<r±{t) =v± = v c ±iv q /2. 



(5) 



The classical component v c is the dot potential gener- 
ated by accumulation of charges with opposite spin. The 
quantum component v„ plays the role of the counting 
field for charge in QD [2l|, |22| • After a number of calcu- 
lations, the CGF is represented by the S'-matrix: 

2 If 
T = .Fq - jjVcVq, Fo = - I du Trln(l - / K) , (6) 



K = l-e ll S\v-)e- ll S{v+) : S- 



Sll &lr 

Srl Srr 
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where 1 is a unit matrix, and A = diag(A, 0) with A= A^ — 
Afl- / = diag(/i, /r) consists of the Fermi distribution 
function f r (u) = l/{exp[(w — fi r )/T] + 1}. When the 
potential v is independent of the magnetic field, the S'- 
matrix is reversible S rr i(v; B) — S r > r (v; —B) 



S rr (v) = 1 — 



iT r + t Ic{ ^T L T R cos (j> - 4f e{v)/2 
A(v) 



Srl(v) = (e*% cf e(v) - ^T L T R ) /A(v), 



ireiV^L^R COS < 



- 1 
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where e(v) = 6D+v — u). The tunnel coupling r = r£ + r/v 
is written with the DOS of the lead as T r = 2nt^.p r . 
Hopping through the reference arm is characterized by 
t rc f = 2'KtLR'JpLpR- It appears that Eq. ([7]) is the 
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CGF for the joint probability distribution of current and 
charge 2jJ,|22|. However, the 'charge counting field' v q 



as a function of A is now determined by coupled saddle- 
point equations: 



v c = {U/2)dT /dv q 



{U/2)dT /dv c 



(8) 



Magnetic field asymmetry in nonlinear transport - The 
saddle-point solution JHJ) captures magnetic field asym- 
metry in the nonlinear transport regime 0, EH- For 



A = 0, Eq. 



possesses a trivial solution: v q 



and 



v c = v* determined by the nonequilibrium Hartree equa- 
tion, v*=U JdwA a {w)/{2iz), 

A a =Y$r + 4 f rv/4)[/ r H - i/2]/|A(«*)l 2 



+ Wv / rTrRsm0[/ L H-/ i? H]/|A(«*)| 2 , (9) 

where f — L/R for r = R/L. Figure Q] (a) shows the 
magnetic field dependence of charge accumulation inside 
the QD, n a = v* /U+l/2. In equilibrium V = 0, n a is an 
even function of the magnetic field. For V ^ 0, because 
the second line of Eq. ([9]) is related to the lack of mirror 
symmetry, the charge accumulation becomes an uneven 
function of AB flux n CT (^>) ^n CT (— <p). 

The average of the charge current is obtained by dif- 
ferentiating the CGF in terms of A. 



dT dJ-Q / dJ-Q dv a 



dv n 



dF 



d{iX) d(iX) J^\dv a d(i\) Ud{iX)J d(i\) ' 

(10) 

where a~c/q for a = q/c. All contributions except T$ 
cancel because of the condition JHJ). Then, the Landauer 
formula with the transmission probablity T = \Slr(v*)\ 2 
is obtained; 



dT a 



A=0 



d(iX) 



T L T R + t 2 ci e(v*) 2 



~[duT{uj) \f L (u) - /r(w)] 



2t Icl e(w*) % /r L r fi cos( 



\A(v*)\ 2 



In the absence of interaction [7 = 0, and thus v = 0, 
the transmission probability is symmetric in the mag- 
netic field [l6j . For finite {/ and V", because of the ac- 
cumulation of magnetic field-dependent nonequilibrium 
charges, the reversibility of S'-matrix breaks down as 
S lr (v*(B);B)^S R l(v*(-B);-B), leading to magnetic 
field asymmetry in nonlinear conductance. 

If we substitute v q = and v c = v* in Eqs. © and 
d?]) , our CGF may be compatible with that in Ref. [7[ at 
the formal level. However, for our case, generally both v c 
and v q satisfying Eq. (jSJ) depend on A. Then the CGF ([6]) 
with Eq. (JSj> , fulfills the FT ([1]) , since if we consider v± 
as variables, .Fo(A, v±; B) = J-"a(— A + iA, v±\ — B) is sat- 
isfied for any v± . Figure Q] (b) demonstrates the FT Q , 
though probability distributions for positive and negative 



magnetic fields are different [inset of Fig. Q] (b)]. There- 
fore, the magnetic field asymmetry does not necessarily 
contradict the FT. 

Nonequilibrium noise.- In the presence of interaction, 
current and charge fluctuations couple in a nontrivial 
manner, which means that we must account for v q care- 
fully. Let us consider the derivative of Eq. @ with re- 
spect to the counting field: 

dv a U d 2 T , U v d 2 T Q dv a , 



d(iX) 2 dv a d(iX) 2 ( jr' cq ® v a 9v a ' d(iX) 



E 



Ugg' 9 2 Tq 

-f 2 dv a > d(iX) 



(12) 



For A = 0, which implies that v q = and v c = v*, 
the four components are U cc \x=a = USnn, U cq \\=o = 
U C q\\=o — U and U qq \\ = o = 0. Coulomb interaction is 
screened U = U/(l — Uxnn), because the right-hand 
side of (|lip contains the derivative of v a itself. The bare 
density-density response function \nn and the density- 
density correlation function (charge noise) Snn are given 
by Xnn = d 2 T /dv c dv q \x =0 /2 = dn a /de D and S NN = 
d 2 T ldv 2 q \ x={) /2. Then, with the help of Eq. {HJ, the 
derivative of Eq. (JTUJ) and the full form of the nonequi- 
librium current noise read as follows: 



d 2 T 



d 2 Fo 



' i \ ' / \ i ai d(iX)dv a 2 d(iX) dv a > ' ' 

((I 2 )) = 2 (S n + 2S IN U xin + Xin 2 S nn U 2 ) , (14) 

where the bare current-density response and the current- 
density correlation functions are xin = (9 ((I)) /den) /2 
and Sin — d 2 To/d(iX)dv q \\—o/2. The current-current 
correlation. 



Sn 



1 d 2 f 



2d(iX) 2 



= i fcLuT(Lu)[f L (Lo) 

A=0 2n J 



-2f L (u)f a (u)] - T(u,) 2 [/lM " fR^)f 



is given by the quantum noise formula for noninteracting 
systems [23[ with the self-consistent potential v*. The 
second and third terms of Eq. (fH)) are the result of in- 
teraction out of equilibrium, since in equilibrium, the av- 
erage current vanishes and consequently, XiN = 0- In the 
absence of the reference arm t re f = 0, Eq. (fT4|) reproduces 
the theory of the noise for the nonequilibrium Anderson 
model in the Hartree-level approximation [2^ |. For this, 
the counting filed of QD charge v q is crucial. 

Nonlinear transport coefficients.- Now we come to the 
relations among the third-order transport coefficients ((3|) . 
First, the bare parts vanish: di\~ n d%!Fa(0, B)\a=o = 
(n = 0, 1, 2, 3). Then, the skewness, following the deriva- 
tive of Eq. (fTB")) in terms of A reads 



L 3 = 6 U^S* x7in > Xii, n = d S n /de D , (15) 
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where xii,n is the linear response of the noise. The su- 
perscript eq. specifies that A is fixed at 0. Equation (fl5|> 
reveals that equilibrium skewness is caused by the inter- 
action. The other transport coefficients are calculated in 
the same manner. 
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where xni — dn^/dA and XijN = 9 Xin /dA. Figured] 
(a) shows the AB flux dependence of third-order nonlin- 
ear transport coefficients. We observe finite skewness for 
4> [panel (a)]. It appears that the coefficients be- 
have independently. However, as shown in panel (b), an 
extension of the Onsager relation (|3|) is satisfied perfectly. 
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FIG. 2: (a) Aharonov-Bohm flux dependent third-order non- 
linear transport coefficients and (b) the extension of Onsager's 
theorem. The parameters are the same as those in Fig. [T] 

We note that our results can be obtained using the 
Hartree approximation based on the nonequilibrium self- 
consistent ^-derivable approximation 
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In this 

scheme, the Keldysh generating function consists of an 
infinite number of closed diagrams, each of which satisfies 
the symmetry ([T]), as shown in Ref. 0. 

Summary- We studied the full counting statistics of 
a quantum dot Aharonov-Bohm interferometer and have 
developed a novel Hartree approximation, which satis- 
fies the fluctuation theorem and describes magnetic field 
asymmetry in the nonlinear transport. We have also 
shown that equilibrium skewness as well as the asymmet- 
ric component of nonlinear conductance are the result of 
Coulomb interaction. These satisfy the extension of On- 
sager relations ([3]) which may be measured by the 
currently available experiments [12[. 
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